General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



-<*“ / 




5- *> 

>A% 

vXv 


w\ 

\\%X 

<L £ % 




(E83-1C057) MAGSAT A NOM AL K FIELD IMEFSION 
AND INTERPRETATION EC E THE US Final Report 
(Business and Technological Systens, Inc.) 

67 p HC AQ4/MF A 3 I CSCL 08N 


G 3/43 


Uncias 

00057 


A v ^ %v 

^ DEC 1982 <3 

(<=0 _ 


RECEIVED £ 

W « *TI FACILITY 





BTS33-82-79/rb 

1001 


\* V *^*' 

, •“ — ♦ Business, iND Technological Systems, Inc 


December 1, 1982 


Business and Technological Systems, Inc. 
Aerospace Building, Suite 440 
10210 Greenbelt Road 
Seabrook, Maryland 20706 


FINAL REPORT 

MAGSAT ANOMALY FIELD INVERSION 
AND INTERPRETATION FOR THE U.S. 


by 

M. A. May hew 


for 


National Aeronautics and Space Administration 
Goddard Space Flight Center 
Greenbelt, Maryland 20771 


Business and Technologigh Systems, Inc 


ABSTRACT 


Long wavelength anomalies in the total magnetic field measured by 
Magsat over the United States and adjacent areas are inverted to an equiva- 
lent layer crustal magnetization distribution. The model is based on an 
equal area dipole grid at the Earth's surface. Model resolution, defined 
as the closest dipole spacing giving a solution having physical signifi- 
cance, is about 220 km for Magsat data in the elevation range 300-500 km. 
The magnetization contours correlate well with large-scale tectonic prov- 
inces. A higher resolution (200 km) model based on relatively noise-free 
synthetic "pseudodata" is also presented. 

Magnetic anomaly component data measured by Magsat is compared with 
synthetic anomaly component fields arising from an equivalent source dipole 
array at the Earth's surface generated from total field anomaly data 
alone. It is found that the synthetic components fit the component data 
regardless of the dipole orientation assigned to the equivalent sources and 
of the dipole spacing. Tentative conclusions are 1) over the U.S., vector 
anomaly fields can be determined to the accuracy of the measurements from 
the total field anomaly data alone, and 2) the equivalent source technique 
is not useful for determining the direction of large-scale crustal magnet- 
ization. 

An excellent inverse correlation between apparent magnetization and 
heat flow in the western U.S. is demonstrated. A new regional heat flow 
map is presented and compared with published maps. A preliminary model of 
magnetic crustal thickness and regional heat flow based on the magnetiza- 
tion model is also presented. The heat flow maps derived indirectly from 
Magsat data "see" nearly all the important thermal anomalies evidenced in 
previous published maps. Notably, the map predicts high heat flow in 
Nebraska and the Dakotas, suggesting the presence of a "blind" geothermal 
area of regional extent. 
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1.0 INTRODUCTION 


This report summarizes the results of experiments in reduction and 
geologic 'interpretation of Magsat anomaly data for the United States. 
Methodologies previously developed and applied to Pogo data (Mayhew et al , 
1980; Mayhew, 1982a) were applied to the Magsat data with the expectation 
that the lower elevation of the latter would lead to higher resolution 
anomaly maps and magnetization source models. Additionally, investigation 
was made of the use of Magsat vector component anomaly data. 

The analysis proceeds in three phases: 1) data processing, 2) anomaly 
inversion to source models, and 3) physical interpretation. The central 
goal of the analysis is to find a magnetization distribution in a thin 
equivalent layer at the Earth's surface having maximum detail while retain- 
ing physical significance, and giving rise to a synthetic anomaly field 
which makes a best fit to the observed field in a least squares sense. 
Multiplying the magnetization distribution by the laye r " thickness gives the 
distribution of the vertical integral of magnetization in the magnetic 
crust to within an indeterminant ambiguity in level (Mayhew, 1982b). Such 
a model can be transformed to a more physically meaningful model given some 
independent constraints, as described below and in Mayhew (1982b). The 
apparent magnetization contrast (Am) in the equivalent layer is approxi- 
mated using an array of dipoles distributed in equal area at the Earth's 
surface. In most applications the dipoles are pointed in the direction of 
the main magnetic field. This carries the implicit assumption that crustal 
magnetization is dominantly induced or viscous. There are good reasons why 
this assumption should generally be valid (Wasilewski and Padovani , 1980; 
Mayhew, 1982a), but an experiment was carried out, and is described in 
Section 3, involving, the hypothesis that large-scale remanent magnetization 
is of importance in some regions. 

The first phase of analysis, data processing, involves 1) selection of 
a quiet data set based on Kp index, 2) subtraction of a main field model 
to produce an anomaly data set, 3) fitting and subtracting a quadratic 
function from each profile to model very long wavelength external field 
effects in mid-latitudes, and 4) application of a running three-point 
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weighted filter to the profiles to minimize very high frequency 
"instrumental" noise in the data. Figure 1 is a track chart of profiles 
used in the present study; coverage is dense. The vector anomaly data was 
treated similarly to the scalar data. However, it is noted that of the two 
vector data sets available on the Investigator-B tapes, the first being a 
selection of every 80^ point from the full data set, the second being an 
80-point average, the latter was so very noisy as to be useless. 

The least squares mathematics involved in inversion of the anomaly 
data is described briefly in the references given above and in some detail 
in Mayhew and Estes (1980a, b). A key element of this phase is the deter- 
mination of the closest possible dipole spacing giving a "stable" inversion 
to a solution having physical significance. This is accomplished by plot- 
ting the standard deviation of the solution parameters (i.e. the magnetiza- 
tion values associated with the dipoles) against their spatial separation 
for a series of solutions. At separations closer than some critical one, 
the magnetization values become large (positive and negative) and no longer 
contour systematically. Figure 2 shows such a graph for Pogo data, imply- 
ing a resolution limit at about 300 km separation. This is to be compared 
with a similar graph based on Magsat data given later in this report which 
implies a resolution limit approaching 200 km for this lower data. Mayhew 
(1982a) shows how to find a magnetization distribution on a grid twice as 
fine as that at the resolution limit, and this technique has. been used in 
the present work. Again, for comparison with Magsat results, Figures 3 and 
4 show a Am distribution derived from Pogo and the associated magnetic 
anomaly field at 450 km. It will be seen that the Magsat results show con- 
siderably more detail in both kinds of maps. 

For the interpretive phase, one of the most strai ghtforward approaches 
Is to attempt to convert the equivalent layer Am models to models of thick- 
ness variation in a layer of constant magnetization. This is applicable, 
for example, to the case in which magnetic field variations reflect undula- 
tions of the Curie isotherm within the crust (this is our working hypo- 
thesis for the western U.S., as described in Section 4). The methodology 
is described in Mayhew (1982b), and is reiterated briefly here. The verti- 
cal integral of magnetization implied by the Am model is H(M+Am), where H 
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is the equivalent layer thickness (arbitrarily taken to be 40km), and M 
is the level ambiguity in the solution, assumed to be constant. For the 
case of thickness variation h in a layer of constant magnetization p , 
the vertical integral of magnetization is hp . These two cases are indis- 
tinguishable in the anomaly field at satellite elevation, so equating the 
two gives a basic equation 


H(M+Am) = ph. (1) 

For the Curie depth problem, if its depth can be estimated at at least two 
places (by use of thermal models or spectral estimates) the unknown M and 
p can be estimated, and h determined for the whole region from the Am 
distribution. The implied Curie depth configuration can then be used to 
constrain a regional crustal geothermal model. Results of this type are 
given in Section 4. 

To summarize the principal elements of the statement of work for this 
project, they are to develop preliminary and final magnetization models for 
the whole of the U.S., use these to constrain thermal models for the west- 
ern U.S., and investigate the implications of such models’ in a resource 
context. 
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2.0 MAPS 


Figure 5 shows a "trade-off" graph for Magsat data analogous to that 
of Figure 2 for Pogo data. The plot of magnetization parameter standard 
deviation against dipole spacing suggests a resolution limit in the vicin- 
ity of 200 km. Note that around this spacing the standard deviation of the 
fit of synthetic and observed (scalar) fields quickly approaches diminish- 
ing returns at a fit around .75 nT. Figure 6 shows the synthetic anomaly 
in the total field at 320 km due to a 136 km dipole spacing. For this 
spacing the fit of model field and data is about as good as can be achieved 
by the equivalent layer approach; it will subsequently be shown that the 
difference is almost entirely due to very high frequency measurement noise. 

2.1 Interim Magnetization Maps 

Our original magnetization map was made at a dipole spacing of 200 km, 
and was presented at an AGU Magsat symposium by Mayhew (1980). The large- 
scale features of this map can be seen in the earlier maps derived from 
Pogo data, but much more detail is present; some of the detail is clearly 
spurious, however, and it was decided- that the resolution limit had been 
pushed too far. We subsequently backed off to 220 km, and obtained a much 
better map (Figure 7) which we believe represents a valid representation of 
magnetization variation with the maximum detail achievable with the scalar 
data. Good correlations with regional geology are present, as noted in the 
caption of Figure 7. Figure 8 is a free-air gravity anomaly map made by 
averaging data on a 300 km grid, and it is noted that good inverse correla- 
tions exist between the magnetization and gravity anomaly maps over much of 
the U.S. 

2.2 Final Magnetization Map 

Resolution of magnetization maps derived from anomaly data is limited 
largely by noise in the data. On the other hand, the field itself can be 
fit extremely well, to within the high-frequency noise level, by an equiva- 
lent layer model with closely-spaced dipoles. This suggests that substitu- 
tion of relatively noise-free synthetic data from such an equivalent layer 
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model ("pseudodata") for real data might permit inversions to higher reso- 
lution magnetization models. Accordingly, a pseudodata set was generated 
using a 136 km spacing dipole array. This dcta set was then inverted to a 
magnetization model involving 200 km dipole arrays. The result was a good- 
looking map very similar to that of Figure 7 3 with a bit more detail, but 
lacking the very questionable contour features of our original 200 km map. 
A second map was made with dipole locations offset 50 km right and up rela- 
tive to the first map. The two maps are nearly identical, differing only 
in very minor detail, which gives one the impression that the whole proce- 
dure provides a valid sampling of a continuous magnetization distribution. 
Finally, the two magnetization grids were averaged on an intermediate grid, 
giving a nice result: the common features of the two maps are retained, 
while the very minor differences are simply smoothed. The resulting map is 
given in Figure 9. Numerous correlations with large-scale geologic fea- 
tures are present. In the western United States, the pattern corresponds 
to large-scale heat flow provinces (e.g. Lachenbruch ana Sass, 1977), and 
may reflect Curie isotherm undulations. In particular, the Basin and Range 
and R 4 ' Srande Rift are regions of high heat flow, and are delineated by 
magnetization lows (possibly indicating a shallow Curie isotherm), while 
the relatively lower heat flow provinces of the Sierra Nevada and Colorado 
Plateau correspond to magnetization highs. The boundary between the 
Appal achian-Ouachit- belt and the Precambrian craton is marked by a strong 
gradient along its length. The Wichita Uplift, flanked on the north by the 
Anadarko Basin, is marked by a magnetization gradient. It is interesting 
that in the western midcontinent the boundary between Mesozoic/Cenozoic 
cover and older rocks marks a distinct change in orientation of the magnet- 
ization anomalies. The Anadarko, Denver, Williston, and Michigan Basins 
all appear to be associated with magnetization positives. Table 1 is a 
list of the parameters of the final Am model. 
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3,0 EXPERIMENT WITH VECTOR DATA 


While in principle vector component data can be used as input for 
inversions to equivalent layer models in which dipole directions are either 
fixed or solved for, in practice the present component data is too noisy 
for this purpose. However, an experiment was carried out with the aim of 
seeing whether the vector data could be used to detect regions of large- 
scale crustal remanent magnetization. The results are reported in Gall i her 
and Mayhew (1982), and are reviewed in this section. 

3.1 Data Treatment 


After "fine-attitude" processing to mathematically redefine the compo- 
nents in Earth, rather than spacecraft, coordinates, a 13th degree and 
order field model containing time terms was removed from the Magsat vector 
component data. Vector data was used in making the field model. The dif- 
ferences between the field model components and the fine-attitude-corrected 
components are the anomaly components. The vector sum of the anomaly com- 
ponents we term "delta-F". The anomaly in the total field ("delta-B") is a 
scalar quantity given by the difference between the magnitude of the vector 
sum of the measured components of the total field and the magnitude of the 
field model vector. Delta-B defined this way is very nearly identical with 
the projection of delta-F in the total field direction. While delta-F is 
uniquely .determined by the anomaly field components, delta-B is not. The 
delta-B data is corrected before input to the inversion program by fitting 
and subtracting a quadratic function in the manner described by Mayhew 
(1982a). The quadratic simulates the very long wavelength fields due to 
magnetospheric ring currents, which are invariably present to some degree. 

3.2 Induced vs. Remanent Crustal Magnetization 

Orienting the dipole sources along the main field direction carries 
the implicit assumption that magnetization by induction (or by viscous 
build-up) dominates over large-scale remanent magnetization in the crust as 
a whole. While several compelling arguments can be made for why this 
should be so, such has not been demonstrated. In applying the equivalent 
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source technique, the dipoles can be oriented in arbitrary directions, and 
a good fit to the observed delta-B field can still be achieved. Thus, the 
delta-B field alone provides no information about the direction of magneti- 
zation in the magnetic source layer. There has been considerable specula- 
tion, however, that use of both scalar and vector anomaly data might pro- 
vide such information. We thought that if the assumption that crustal 
magnetization is dominantly in the main field direction is valid and we 
inverted delta-B data to a magnetization solution, with sources at the 
critical spacing (220 km) and oriented in the main field direction, that 
the associated synthetic vector field should fit the measure vector anomaly 
components. We found this in fact to be the case: the equivalent source 
vector components fit the measured components along arbitrarily selected 
satellite profiles within the noise in the data (Figures 10 and 11), 
provided an appropriate quadratic function is removed from each component 
separately. We initially took this to be support for the induction 
hypotheses. We then tried two experiments in which the dipoles were 
oriented in odd directions. In the first, all were pointed in declination 
45°E, inclination 35°; in the second, all were pointed with declination 
45°E, inclination -35°. Again delta-B was inverted to a magnetization 
distribution. We expected to be able to come back with a fit to delta-B 
but not the components, but in fact in both cases the fit to the components 
was very nearly the same as the original fit. As a corollary experiment, 
we inverted delta-B to a magnetization distribution for a much smaller 
dipole spacing (136 km, an arbitrary value). Sources were pointed in the 
main field direction. Because the magnetization values are large and 
irregular at this spacing, we expected to fit delta-B slightly better than 
above the critical spacing, but not the components. Again, however, the 
fit was comparable with previous experiments, although in fact the fit of 
delta-B is very slightly better at this spacing than at the resolution 
limit. All experiments were repeated for several groups of profiles over 
three different areas of the U.S. with the same result. Figures 10 and 11 
show two examples of the fits. In each column of these figures is shown, 
from bottom to top, the z, y and x components, the anomaly in the total 
field ("delta-B"), and the magnitude of the' anomaly vector ("delta-F"). 
Small dots are raw anomaly data, while the smooth solid line is a quadratic 
function fit to this data; the scale for these plots is on the right. The 
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residual between them is shown as an asterisk. The quadratic is that which 
brings the residual into best agreement with the equivalent source synthet- 
ic field (irregular solid line); the scale for these is at the left. 
Dipole sources are not placed within about 32° of the magnetic north pole 
to avoid auroral effects, so that the equivalent source fields tend to zero 
at the higher latitude. 

The rather surprising result described above suggests that, with the 
equivalent source technique, the vector anomaly field is uniquely determin- 
ed by delta-B, and that it cannot be used to detect large-scale remanence. 
We hope at a future time to theoretically verify this empirical result, 

3.3 C omments on Data Characteristics 

The equivalent source field provides a reference against which the 
various sources of "noise" in the measured data can be evaluated (Figures 
10 and 11). Clearly, the component data is much noisier than the delta-B 
data, although it is well within the mission specifications. Noise is of 
several kinds. First is the very long wavelength ring current effect 
modeled as a quadratic function. The 'three quadratic functions associated 
with the three components define a vector which may be of interest in 

future studies of the ring current field. Second, some of the profiles 

• 

show a strong perturbation at high latitude approaching the auroral zone. 
These higher latitudes have been avoided in determining the equivalent 
source solution, so that the synthetic field goes to zero there. Third, 
the component data commonly shows local divergence from the synthetic field 
followed by little jumps. This is due to the accuracy limitations on the 
fine altitude determination. It is suggested that the equivalent source 
technique could help improve the accuracy of the attitude determination. 
Fourth, the component data commonly shows small undulations about the syn- 
thetic field which are not consistent from profile to profile and undoubt- 
edly are small-scale external field effects. 
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4.0 INTERPRETATION 


The interpretive phase of the present study was mainly an investiga- 
tion of the relation between 4m anomalies and heat flow for the western 
United States. Mayhew (1982b) showed that an inverse relation exists for 
regionally averaged heat flow and low-resolution Am variations inferred 
from Pogo data. In the first part of this section it is shown that a simi- 
lar relation exists for Magsat data on a finer scale. The second part 
describes an investigation of'whether Am maps can be used in a more quan- 
titative fashion to accurately predict heat flow and Curie depth. 

4.1 Heat Flow Map 


The production of regional heat flow maps is problematical for a num- 
ber of reasons. Heat flow data sets are inherently "noisy" because of 
various transient, ground water, topographic, and other effects. Measure- 
ments are irregularly distributed and of variable quality. Conductive heat 
flow is difficult to separate from convective effects due to local move- 
ments of hydrothermal fluids or magmas in active areas. 

Curie isotherm undulations should be reflected in regional heat flow 
variations. To the extent that Am anomalies are due to Curie depth vari- 
ations, -they can be used along with heat flow measurements to guide the 
construction of regional heat flow maps. 

An up-to-date heat flow data set was used to make two equal -area aver- 
age grids slightly offset from one another. Grid spacing is 200 km. The 
averaging smoothes local anomalies, and the two offset grids gives a kind 
of running average presentation. The grid average values were printed 
together on a map and contoured using the Am anomaly map (Figure 9) as a 
guide. The result is given in Figure 12. The idea here is that if the 
above hypothesis is correct it should be possible to contour the average 

* 

heat flow data in such a way as to be inversely related to Am. Figure 12 
shows that this is indeed possible for nearly the whole of the western 
U.S., the only problem area being parts of the Basin and Range. The diffi- 
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culty in this region is likely to be convective effects, hydrothermal in 
the northern areas, hydrologic in southern Nevada, as discussed by Mayhew 
(1982b); note that for these problem areas adjacent average heat flow 
values are inconsistent with each other, probably because the values are 
biased by local extreme measurements due to convective phenomena. Notable 
features of the heat flow map of Figure 12 are as follows. 

The principal heat flow lows are the Sierra Nevada-Baja California 
belt and the zone extending from the Colorado Plateau north through western 
Wyoming. The strongest heat flow high is in the Yellowstone area, broaden- 
ing south to include the whole of the Basin and Range and its extension 
south of the Colorado Plateau. The "Battle Mountain high" of northern 
Nevada, a conspicuous feature of all heat flow maps based on measurements 
is reflected in the contours, but the "Eureka low" of southern Nevada is 
not; the latter has been attributed to regional ground water movement 
(Lachenbruch and Sass, 1977). Likewise* the Am contours suggest much 
higher regional conductive heat flow along the eastern margin of the Basin 
-VKi ttange than is indicated by measurements. The Yellowstone high is 
inferred to extend northwest in eastern Idaho and Washington; the heat flow 
data support this inference. 

A significant heat flow high is inferred for Nebraska and South 
Dakota, possibly extending through North Dakota. While some high heat flow 
values of the area support this idea, there is no obvious surface expres- 
sion of tectonic or volcanic activity usually associated with regions of 
high heat flow, except for the Tertiary volcanics of the Black Hills up- 
lift, suggesting the possible existence of a "blind" geothermal area of 
regional extent. Hot ground water is known from some parts of the region. 

A heat flow high associated with the Rio Grande Rift is strongly ex- 
pressed in the Am contours as in previous maps derived from Pogo data 
(Mayhew, 1982b). The Am map suggests that the heat flow anomaly extends 
well into southern Wyoming, but this is not supported by measurements. The 
Am contours also suggest higher heat flow in southwestern-most New Mexico 
than is perhaps suggested by measurements. 
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The Inference by Mayhew (1982b) that the thermal anomaly associated 
with the rift turns to a southeast trend in northern Mexico is not support- 
ed by the Ain contours or the average heat flow values. It is curious 
that the Am contours in northern Mexico have a southwest trend quite 
oblique to regional tectonic trends, but seemingly consistent with the heat 
flow pattern. 

Figures 13a and 13b are two recently published heat flow maps. Figure 
13a is taken from Blackwell (1979), and is based on measurements; 13b from 
Swanberg et al (in press) is based principally on silica geothermometry. 
There is a remarkable correspondence between features of these maps and the 
heat flow map presented here. A comparison is tabulated in Table 2. 

Figure 14 shows the distribution of major mineral localities in the 
western U.S. with respect to the Am distribution. One gross generaliza- 
tion can be made, that with the exception of the Sierran localities, nearly 
all of the mineral deposits lie within or along the margins of Am lows. 
This is to be expected if the lows are associated with elevated heat flow 
which is in turn related to young magmatic or hydrothermal activity. It is 
suggested that more detailed investigation of these relationships in local 
areas might profitably be undertaken. 

4.2 Heat Flow/Curie Depth Models 

In this section is described a quantitative approach to modeling heat 
flow and Curie depth variations using the Am model derived from Magsat 
data, under the hypothesis that Am variations reflect variations in 
depth-to-magnetic bottom. The basic equation is (1), from which variations 
in magnetic crustal thickness h are computed from variation in Am. 
Estimates are needed for the parameters M and y in equation (1). These 
are gotten by fitting a line to a plot of estimates of range of h derived 
from spectral depth estimates, depth to 550°C (the assumed Curie tempera- 
ture) given by thermal models constrained by surface heat flow measure- 
ments, or, in the case of the Sierra Nevada, seismic estimates of depth to 
mantle, against range of Am for the corresponding areas. The Moho is 
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assumed to be the magnetic bottom for the Sierras, since the very low heat 
flow values measured along the axis of the belt imply depths to 550° well 
into the mantle, which is assumed to be non-magnetic (Wasilewski et al , 
1979). The plot is given as Figure 15. Sources of the values plotted are 
listed- in Table 3. The ranges of values are rather large, but give an 
overall trend in the right direction. The dashed line shown implies values 
for M and y of 0.65 and 1.0 A/m, respectively. These values can be 
used to compute h for the whole map area, and this is shown for the above 
model by the contours in Figure 16. Also shown for comparison are esti- 
mates of depth to 550°C. given by a one-dimensional thermal model using a 
set of average heat flow values. The model assumes constant surface heat 
production decaying exponentially with depth, and constant thermal conduct- 
ivity. The equation is given by Mayhew (1982b). The agreement is fair of 
some areas, poor in other, notably much of the Rio Grande Rift. Assuming 
the model depth values contoured in Figure 16 represent an isothermal sur- 
face, the same thermal model can be used to compute predicted heat flow. 
This is done in Figure 17 (contoured values) and compared with the grid 
average measurements (printed values). Again, while favorable agreement 
exists in some areas, there is disagreement in others. The exercise 
undoubtedly suffers from problems arising from many simplistic assumptions 
on both sides, data limitations, and blanket treatment of a very large 
region. In fact it is somewhat remarkable, and even encouraging, that 

there is any agreement at all . It is worth listing the assumptions inher- 
ent in the analysis. 

• Magnetization is constant throughout the magnetic crust. 

• The Am model level ambiguity is constant throughout the map area. 

• Am variations reflect only magnetic crustal thickness variations, 
which are completely resolved. 

• Curie temperature is constant. 

• Heat production in the crust decays exponentially with depth from a 
constant surface value. 

• Thermal conductivity is constant throughout the crust. 

• Gridded heat flow values accurately represent area averages. 

• Heat flow is by conduction only. 

• One-dimensional heat flow models give a valid first-order approxi- 
mation. 
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While one can have no great confidence in the quantitative estimates 
of regional heat flow given by Figure 17, it is suggested that the method- 
ology is promising provided that 1) more limited areas are treated, 2) 
better estimates of model parameters are obtained, and 3) a critical evalu- 
ation of available heat flow measurements and of the importance of the 
convective component is made. Geothermal resource potential maps can be 
produced from a good thermal model derived from a Am model. For example, 
the theoretical conductive heat content per unit area in the outer 3 km of 
the Earth implied by models such as those of Figures 16 and 17 is given by 
the following equation (Diment et al , 1975) 


TJ - C ll Tdz, 


where TJ is heat content in a 3 km column of unit area, C is volumetric 
specific heat (assumed constant and equal to 0.6 cal /cm 3 °C), T is tem- 
perature, z is depth. 
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FIGURE CAPTIONS 


Fig. 1 Magsat satellite tracks indicate data coverage for this study. 


Fig. 2 Trade-off between (1) fit of equivalent source magnetic anomaly 
field arising from equal area source array at Earth's surface and observed 
field at satellite elevation and (2) 'stability' of inversion as measured 
simply by standard deviation of solution parameters. Optimal solution 
inferred to be at source spacing of about 300 km. Graphs for eastern and 
western halves of map area shown separately; sd means standard deviation. 
From Mayhew (1982a). 


Fig. 3 Apparent magnetization contrast in 40 km thick layer obtained from 
inversion of satellite magnetic anomaly data. Model is series of 300 km 
dipole grids computed separately but machine contoured together without 
smoothing on 150 km grid. Contour interval is 0.1 A/m. For simple model 
of uniformly magnetized magnetic crustal layer, magnetization values are 
proportional to magnetic layer thickness. From Mayhew (1982a). 


Fig. 4 Equivalent source representation of magnetic anomaly field at 450 
km elevation. Equal area source distribution at 200 km spacing. Contour • 
interval 0.5 nT. From Mayhew (1982a). 


Fig. 5 Trade-off as a function of dipole spacing (kilometers) between 1) 
the fit of the equivalent source magnetic anomaly field to the field 
observed at satellite elevations (tenths of nT) and 2) 'stability* of 
inversion as indicated by standard deviation of solution parameters (tenths 
of A/m). Optimal solution inferred to be source spacing of about 220 km. 
SD means standard deviation. From Mayhew and Galliher (1982). 


Fig. 6 Equivalent source representation of magnetic anomaly field at a 
height of 320 km. Source spacing is 136 km. Contour interval InT. Albers 
equal area projection. From Mayhew and Galliher (1982). 


Fig. 7 Apparent magnetization contrast in a 40 km thick layer. Distribu- 
tion is obtained by inversion of Magsat total field anomaly data. Model is 
series of staggered 222 km - spacing dipole grids computed separately, but 
machine contoured together without smoothing on 111 km grid. Contour 
interval is 0.1 A/m. Dashed lines indicate generalized tectonic province 
boundaries which are geographically related to the magnetization distri- 
bution. Letters represent the following; S: Sierra Nevada block, BR: 
Basin and Range province, CP: Colorado Plateau, Rr Rio Grande Rift, SR: 
Snake River Plain, P: western boundary of Great Plains, A: Appalachian 
Ouachita front, W: Wichita uplift, K: Kentucky anomaly, LS: lake 
Superior Syncline. Boundary between Mesozoic/Cenozoic cover and older 
rocks is indicated by short dashed line. Albers equal area projection. 
From Mayhew and Galliher (1982). 
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Fig. 8 300 km grid average free-air gravity. Contour interval 10 mgal. 


Fig. 9 Apparent magnetization contrast (Am) in 40 km layer derived from 
pseudodata set as described in text. Contour interval for 9a is 0.05 A/m. 
Contour interval for 9b is 0.2 A/m. Some geologic provinces labeled as in 
Figure 7. Other features are MCR: Mid-Continent-Rift, Ad: Adirondacks, 
Ab: Anadarko Basin, Db: Denver Basin, Wb: Williston Basin. 


Fig. 10 Comparison of data with equivalent source magnetic anomaly field 
computed from 136 km grid for nearly north-south Magsat profiles over 
eastern U.S. at about 276°E. Horizontal scale is degrees latitude; 
vertical scale is nT. From Gall i her and Mayhew (1982). 


Fig. 11 As in Fig. 10 for two nearly north-south profiles over western 
U.S. midcontinent at about 260°E. From Gal li her and Mayhew (1982). 


Fig. 12 Contours of heat flow averaged over two staggered 200 km grids 
(one grid shown) and guided by Am distribution of Figure 9. Contour 
interval 0.5 HFU. 


Fig. 13 Two published heat flow maps to be compared with that of Figure 
12. Figure 13a from Blackwell (1979), 13b from Swanberg et al (in press). 


Fig. 14 Major western mineral localities from Ander (1981) compared with 
Am distribution of Figure 9. 


Fig. 15 Estimates of magnetic crustal thickness based on aeromagnetic 
spectra, depth to 550°C (assumed Curie temperature) , or, for Sierran 
region, seismic crustal thickness. Dashed line implies true (constant) 
magnetic crustal magnetization of 1 A/m. Sources of estimates given in 
Table 2. 


Fig. 16 Estimates of depth to 550°C. based on one-dimensional thermal 
model and Am distribution shown as contours. Contour interval 10 km. 
Depth to 550° based on same thermal model and observed heat flow averages 
printed. 


Fig. 17 Theoretical heat flow implied by Am distribution and thermal 
model shown as contours. Contour interval 0.5 HFU. Average observed heat 
flow values printed. 
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Table 1 


Parameters for x-y grid (60 x 32 array of dipoles) which over- 
lays 1:7‘5 million scale Albers projection. I and J are the x and y 
positions, respectively, of individual dipoles. LAT and LON are lati- 
tude and longitude (east) . X and Y are the x and y map coordinates in 
inches. DELTA-M is the apparent magnetization contrast value associated 
with the dipoles (tenths of A/m). Flag value 999.00000 indicates un- 
defined value at that location. 
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POSITIVE OR 


CLEARLY "SEEN"? 

LOCATION 

NEGATIVE 

Blackwell 

Swanberg 

this study 

Sierra Nevada-Baja California 

- 

yes 

yes 

yes 

Colorado Plateau-Wyoming Basin 

yes 

yes 

yes 

Basin and Range-Yell owstone 

+ 

yes 

yes 

yes 

Rio Grande Rift 

+ 

yes 

yes 

yes 

Central Montana 

- 

yes 

yes 

yes 

Battle Mountain anomaly 

+ 

yks 

(no) 

(yes) 

Eureka anomaly 

- 

yes 

no 

no 

Cascades anomaly 

+ 

yes 

no 

no 

Northern Great Plains 

+ 

(yes) 

yes 

yes 

Texas Panhandle 

+ 

- 

yes 

no 


Table 2. Comparison of features "seen" or not "seen" in regional heat flow 
maps of Blackwell (1979), Swanberg et al (in press), and this study. 
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Table 3. Source of information on 
magnetic crustal thickness plotted 
in Figure 15. 


1 Sierra (38°N) 

2 Sierra (whole) 

3 Uinta Basin 

4 Utah High Plateaus 

5 Yellowstone 

6 Colorado Plateau 


Eaton (1963) 

Pakiser and Brune (1980) 
Shuey et al (3977) 

Shuey et al (1977) 
Bhattacharyya and Leu (1975) 
McGetchin and Silver (1972) 










